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Abstract 

We give definitions for different types of moving spatially localized objects 
in discrete nonlinear lattices. We derive general analytical relations connect- 
ing frequency, velocity and localization length of moving discrete breathers 
and kinks in nonlinear one-dimensional lattices. Then we propose numerical 
algorithms to find these solutions. 

I. INTRODUCTION 

The search for moving radiationless spatially localized excitations is an interesting topic 
in the field of nonlinear dynamics of systems of many interacting degrees of freedom. 

Many of the integrable one-dimensional models, both discrete and continuous in space, 
posess moving objects, either breathers or kinks QXJ] . The integrability of these models pro- 
vides one with action-angle variables. Tuning the actions as parameters one can continuously 
go from stationary solutions to moving ones (examples are the nonlinear Schrodinger equa- 
tion [|1], the Toda [Q] and Ablowitz-Ladik lattices ||). This fact allows us to think about the 
whole family of solutions as of a particle-like entity. Still one has to admit that the precise 
connection between integrability and possible existence of moving localized objects is not 
known. The reason is the hidden character of the symmetries which provide integrability. 

Another reason for the existence of moving solutions can be some continuous symme- 
try of the Hamiltonian, e.g. the invariance under Lorentz transformations. The Lorentz 



transformation generates moving objects provided the corresponding stationary object ex- 
ists. The system does not need to be integrable, so moving kinks exist for instance in $ 4 
theory in 1+1 dimensions [|1|]. However stationary breathers appear to be nongeneric for 
space-continuous models already in the case of one spatial dimension ||, The reason 
for that lies in the fact that the phonon frequency spectrum Q q (q is a wave vector) of 
small-amplitude (linearized) vibrations of a continuous system is typically unbounded from 
above. Per definition a stationary breather is a time-periodic spatially localized solution of 
the equations of motion. Expanding the solution in a Fourier series with respect to time 
one has to deal with the Fourier components associated to each multiple kQf, of the breather 
frequency Qf, with integer k. The unboundeness of the phonon spectrum leads to an infi- 
nite number of high frequency resonances, which generally prevents from the appearance of 
space-localized time-periodic solutions |J. An exception is e.g. the sine-Gordon system and 
some isolated perturbations of it. 

Space-discrete models can generically allow for the existence of stationary discrete 
breathers |7|], M (see || for an extensive discussion and [Tt| for a list of references). One 



reason for that is that the discreteness of space produces a lower cutoff in the wavelength 
of small amplitude plane waves, and thus a finite upper bound for the phonon spectrum. 
Then one has the possibility to choose frequencies Qb such that kVt^ ^ Q q for all integer 
k. However discreteness in space implies the loss of say the continuous Lorentz symmetry. 
Thus in general there is no clear way how to generate moving breathers out of stationary 
ones. 

To look for moving breather solutions we need to have some good definition of them. 
We define the simplest type of a moving solution as a solution that repeats itself after the 
time T s shifted by one lattice site. Such a solution is a fixed point of the map RGt s where 
G t is the evolution operator in the phase space of the system, R is the translation operator 
that shifts all indices by 1. A little bit more sophisticated but still simple solutions can be 
obtained by considering fixpoints of the map R n G, i. e. solutions that repeat themselves 
after the time T s shifted by n sites. We assume the lattice spacing to be 1 so the velocity V 



is then just n/T s . 

Depending on the boundary conditions at infinity we speak about moving breathers 
or kinks. A trivial example of a radiationless moving discrete kink can be obtained by 
considering identical billiard balls on the line separated by some distance /. Then kicking 
one ball we get an eternal motion where after the n-th kick the n-th ball transfers all its 
energy to the n+l-th one. The interaction potential of balls can be more or less arbitrary 
the only restriction is that it is short-ranged enough, so only two balls interact at each 
moment of time. Below we will in general discuss boundary conditions corresponding to 
moving breathers, i.e. the lattice is asymptotically in the same groundstate no matter what 
direction from the center is chosen, given only that the distance from the center becomes 
infinitely large. 

In the next section we will outline necessary conditions o fexistence of moving breathers. 
In section III examples are given, and numerical calculations of solutions are presented in 
section IV. In section V we relate our results to the work of others. 

II. NECESSARY CONDITIONS OF EXISTENCE OF MOVING BREATHERS 

Since a stationary breather is characterized by an internal frequency Qf,, we have to incor- 
porate this timescale into the definition of a moving breather. Consider a one-dimensional 
lattice, describing the interaction of degrees of freedom associated to each lattice site. Each 
degree of freedom is given by a pair of canonically conjugated variables (e.g. displacement 
and momentum) labeled with the site index. Call one of those variables u n (t). We define a 
one-frequency discrete moving breather solution as 



Here F(x,y) is a function period 2n periodic with respect to x and localized with respect 
to y: 



u n (t) = F(Q b t,n-Vt) 



(1) 




(2) 
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If T s and 2ir/Qb are commensurate so kT s = I2ir/Qb, where k and I are integers, then such a 
breather repeats itself after time kT s shifted by k sites and belongs to the simplest moving 
breathers defined above. In the general noncommensurate case the breather will never repeat 
itself although coming arbitrary close to it. 

In the same manner breathers having two or iV internal frequencies can be defined. This 
hierarchy of objects incorporates everything that we intuitively percept as an object moving 
through the lattice. 

Thinking of moving breathers in terms of fixed points allows to define other interesting 
objects on a discrete lattice. Consider a fixed point of some general map Gt s X where X is an 
element of the lattice symmetry group. If X is the identical transformation we get stationary 
breathers. The translation operator gives us moving ones. For one-dimensional lattice the 
only symmetry group element left is the reflection, which gives us reflector-breathers, which 
mirror themselves after time T s . Higher dimensional lattices provide more choices namely 
taking a rotation as X we get rotation-breathersQ , and then taking as X a superposition of 
rotation and translation we get "walking" breathers. 

It is clear from the beginning that looking for discrete breathers in terms of exact analyt- 
ical solutions is a hard task. Approximate methods like a rotating wave approximation even 
if justifiable would turn us away from the phase space of the system and thus do not help 
too much to understand what a moving discrete breather is. A productive way is to look 
at moving breathers from a general point of view and although not solving any particular 
problem to exactly find model independent relations that all the moving breathers should 
fulfil. 

We derive these relations below considering the tails of a moving breather where the 
motion can be considered linear because of the small amplitude of oscillations. 

Let us consider a moving breather with one internal frequency as defined above in ([[]). 



1 Not to be confused with rotobreathers, see jjTl|], [12], [13|, [14]. 
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We write the function F(x, y) in a Fourier series with respect to x : 

F(x,y)=J2e lkx f k (y) • (3) 
k 

Inserting this ansatz into the Hamiltonian equations of motion we obtain coupled equations 
for the functions f k (y)- In the spatial tails of the breather these equations decouple with 
respect to the label k. Let us consider some particular value of k with frequency fl k = kQ b . 
The equations are linear so we seek for a solution in exponential form: 

h(y) = e-^ . (4) 

Here q k is a complex number. The u n (t) term corresponding to a given k-th harmonic then 
takes the following form: 

u nk {t) = e**-^ . (5) 

Here z = Q + Vq k . Using the equations of motion we finally obtain 

z = G(q k ) , (6) 

where G(q) = Q q is the dispersion relation of the system analytically continued to the 
complex plane. 

With the definition of z we have 

%z = VZq k , (7) 

where (S^) -1 = A is the localization length of a given harmonic, 9ffc = Q + V^R,q k , ft = kVt b . 

This equation connects the frequency, velocity and the localization length of a given 
harmonic with each other. For the particular case of a breather without internal frequency 
Slf, = (shock wave) we obtain the following relation between the velocity and the localiza- 
tion length: 

Vq k = G(q k ) . (8) 
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This equation shows that for a shock wave with a given velocity only a discrete set of 
complex wave vectors is allowed. For the linear dispersion relation G(q) = Cq only the 
velocity V = C is allowed and then q can be arbitrary. For the parabolic dispersion relation 
G(q) = Aq 2 q is always real namely q = V/ A so no localization is possible because the 
localization length is the inverse of the imaginary part of q. For the quadratic dispersion 
relation with a cutoff at finite frequency 

Vq = G(q) = Aq 2 + D (9) 

we have a complex root when (V) 2 — 4AD < 0. Therefore a cutoff in the dispersion relation at 
finite frequencies in the dispersion relation is a necessary condition to have localization of a 
moving solution. Note that such a cutoff can be caused either by a gap in a space-continuous 
system or simply by a space-discrete system. 

For moving breathers with nonzero internal frequency a necessary condition of existence 
is obviously that all harmonics fk{y) are localized namely for any Q = kQf, there is a solution 
of equations (|7|) having a complex q k . 

Therefore we obtain that for any nonlinear lattice the frequency of the breather, its 
velocity and the localization length of its harmonics are not independent but connected by 
the equation (0). 

III. EXAMPLES 

In this chapter we will consider examples of systems to illustrate the general considera- 
tions from above. 

A. Space-continuous models 

1. Klein-Gordon models in 1+1 dimensions 
Consider the partial differential equation for the field U(x, t) 
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U u = -U + CU xx -F nl {U) , (10) 

where the function F n i(x) if expanded in a Taylor series around x = contains only nonlinear 
terms in x. Examples are e.g. the $ 4 Klein-Gordon equation with F n i(x) = x 3 or the sine- 
Gordon equation with F n i(x) = sin(x) — x. 

Let us search for a moving solution in the form 

U{x,t) = J2 u k(x-Vt)e m . (11) 
k 

Inserting (0) into (0), assuming that the amplitude of the solution is small at large dis- 
tances from some center, skipping the nonlinear terms in (^) we obtain for each integer 
k 

(v 2 - c)^- - 2iVkn b ^ + (i - k 2 n 2 b )u = o . (12) 

az A dz 



For simplicity we skip the k index and define u = u^, = kil^, z = x — Vt. To solve ( |T2"|) 
we make the ansatz u(z) ~ e Az . Note that A = iq^ (see (f|)). Decomposing A into real and 
imaginary parts A = R + il (R, I real) we find 

(V 2 - C) (R 2 - I 2 ) + 2vni + 1 - tt 2 = (13) 
(V 2 - C)2RI - 2VQR = . (14) 

The physically relevant parameters of our solution are the velocity V and the exponent of the 
spatial decay R characterizing the localization length of the object. The physical frequency 
describing the true oscillations is given by uj = Q — VI. Solving (pT3|) we find 
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I 2 = V - 

c 



1 



C-V 2 



-R 2 



(15) 



u = — I. (16) 

It follows that V 2 < C (Lorentz invariance of ( |i~0|) and V 2 > C — 1/R 2 . These two curves 
define the allowed region in the {-R 2 , V 2 } plane of possible solutions, which is shown in Fig.l. 
The velocity of any moving object has to be below the speed of light. For large (weakly 



localized) objects with R 2 < 1/C the solution will always have some nonzero frequency u, 
i.e. there is a gap in the allowed frequency spectrum containing to — 0. For small (strongly 
localized) objects with R 2 > 1/C the gap closes, and one can always design a tail solution 
with zero frequency uj = (this corresponds to the lower bound of the allowed region in 
Fig.l). The case of zero frequency is nothing but a tail (or front) of a shock wave. Note 
that there are no restrictions with respect to the value of R 2 , so the tail solutions can be 
infintely stronlgy localized in space. 

So far discussed the solutions without checking whether the initial ansatz (|TT|) can be 
completely fulfilled in the tails, i.e. for all integers k. Note that in ansatz ([ll]) we parametrize 
the solution using Qf, and V. To answer that question, we can argue in the following way. 
Suppose we choose a point in the {-R 2 , V 2 } diagram in the allowed region. This point gives 
us a set of values for u and /. Suppose that k = 1. Then we obtain some unique value for 
Qb. Now we can go the inverse way and say, that for that value of Qb and V we find the 
corresponding values for R and /. But what now for other values of kl Since increasing k 
we increase u>, we have to check whether at a fixed value of V we can realize any value for 
uj in Fig.l by changing R throughout the allowed region. The answer is no. Indeed fixing 
V we always obtain a finite line segment in Fig.l. On the right end of this segment uj = 0, 
and on the left uj < oo, with no singularities in between. Thus we can never realize ansatz 
(p|). So we conclude that in general moving breathers do not exist in Klein Gordon field 
problems. The only possiblity is to have a shock wave (or kink), i.e. to set = I = uj = 0, 
which is possible. 

It seems that we are in conflict with the well known result, that stationary and moving 
breathers exist in the sine-Gordon PDE. But keep in mind, that first the sine Gordon PDE is 
integrable. Secondly if one takes the stationary breather solution for that case and expands 
it in a Fourier series with respect to time, inserts the sum into the equations of motion and 
looks at the behaviour of the Fourier components in the spatial tails of the solution, it appears 
that their decay is not described by the linearized equations of motion. In this nongeneric 
(since integrable) case some nonlinear terms in the equations of motion have to be kept to 
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explain the final exponential decay of breather harmonics which resonate with the phonon 
band. This corresponds to the following: suppose you solve a differential equation with 
inhomogeneous terms. Then the full solution is a sum of the general homogeneous solution 
and a particular inhomogeneous solution. In our context the homogenenous part comes from 
the linear terms, and the inhomogeneous from some nonlinear terms. The high symmetry 
of the sine Gordon equation miracuously ensures the vanishing of the coefficients in front of 
the (nondecaying) homogeneous part. However as already discussed most perturbations of 
the sine Gordon equation, by destroying those hidden symmetries, destroy the breather. So 
although this is a subtle point, we don't see any contradiction to our generic statement. 



2. (Non)linear Schrddinger equation 

Consider the partial differential equation for the complex field ty(x,t) 

* = i(CV >xx + F(V)) . (17) 

Again we search for a solution in the form = 0(x — Vt)e int . Repeating the same 

procedures as in the previous case we arrive at the equations 

/ = ^, (18) 

C(R 2 -I 2 ) . (19) 



Since ( fPTD is not Lorentz- invariant, we do not find restrictions on the choice of the velocity 
V. In fact we find no restrictions at all, the whole parameter space {R, V} is allowed. 

B. Space-discrete models 

1. Discrete (non)linear Schrddinger equation 

The equations of motion are given by 

* n = i[|* n |^ n + C(* ft _i + *„+!)] . (20) 



Here the nonlinear term is characterized by some \x > 0. Again we search for a moving 
solution in the form \l/ n = <f>(n — Vt)e . Due to the nonlocality of the difference operator 
(as compared to the differential one) the differential equation for (p(z) contains now retarded 
and advanced terms (the tail is again only considered, nonlinear terms are neglected): 

-V4>'(z) = i[-n<f>(z) + c(4>(z + i) + 4>(z-i))} . (21) 

Skipping the intermediate calculations we arrive at the result 

VR 

I = arcsin , (22) 

l2sinhi?J ' v ; 

uj = 2 coshi? cos/ . (23) 

A necessary condition is thus V < 2smh(R) / R and is shown in Fig. 2. In the region of 
possible solutions for each pair of {R, V} we now have two solutions due to the periodicity 
of (|22| ) , (|23|) in I - to each solution with a given value of Ii — I we can construct a second 
solution with I 2 = tt — I. For any given value of R there is a finite upper bound on the 
value of V. However with increasing R the threshold value of the upper bound of V also 
increases. Again we can find infinitely strongly localized and infinitely fast moving tails. 

2. Klein- Gordon chains 

These models describe the dynamics of atoms on a substrate and interacting with each 
other: 

u n = -au n - C(2u n - w n _i - u n+ i) +F n i(u n ) . (24) 
As in the field case we search for a moving solution in the form 

u n (t)=Y J Mn-Vt)e ik ^ t . (25) 

k 

In the tails of the assumed existing solution we obtain 

y2 fMz)_ _ 2ikQb ydA^zl = {n2 _ a _ 2C )A k (z) + C(A k (z + 1) + A k (z - 1)) . (26) 
az z az 
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Repeating the intermediate calculations as above we arrive at the following equations (note 
that we skip the index k, so below to = kQ b — VI): 

C 

to = —-rr-jr sinh/2 sin/ , (27) 
V R 

C 1 2 

sinh 2 i? sin 2 / = V 2 R 2 + a + 2C(1 - coshi? cos/) . (28) 



VR\ 

The allowed region in the parameter space {R, V 2 } is similar to the one of the discrete 
nonlinear Schrddinger case (Fig. 2) but also more complicated. Solutions exist below a certain 
line V 2 (R). This line consists of two parts. For cosh/? > 2(7/ (a + 2C - y/a 2 + 4aC) the 
line is given by 

(a) : V 2 = — (2CcoshR-a-2C) . (29) 
R 

For smaller values of R we have line (b): 

1 sinh 2 /? 



(b) : V 2 = -— — (a + 2C- Va 2 + AaC) . (30) 

In this second case the continuation of line (a) (|29|) separates solutions with multiplicity 4 
(to the left) from solutions with multiplicity 2 (to the right) as shown in Fig. 3. 

Let us fix any value of V. This value defines us some half-infinite line segment of allowed 
solutions in {R, V 2 }. For uj we find < uj < oo. Thus for any V ^ we can generate a 
breather solution in the tails, for any value of fi&! Surprisingly the problem of resonances, 
as in the case of a stationary breather, does not appear on that stage. 

We can also find a solution to f2& = uj = I = 0, i.e. we can again generate shock waves 
(or kinks) with any velocity (in the tails). They correspond to solutions on line (a). 

3. Acoustic chains 

This case is obtained by performing the limit a — > 0. The phonon spectrum is acoustic. 
In that case solutions exist below line (a) (29[), which now extends down to R = (Fig.4). 



All solutions are of multiplicity 2. In contrast to the previous case for velocities V 2 < C we 

11 



can generate moving breather tails only for frequencies above some threshold value. This 
gap value shrinks to zero and remains zero as the velocity is increased above V 2 > C. Also 
shock waves (kinks) which correspond to all points on line (a) can be again generated. Again 
in contrast to the previous case, these shock waves or kinks must have velocities V 2 > C to 
be generated. 

C. Some discussions 

Let us summarize the results obtained so far. We assumed the existence of some moving 
localized object (breather) which is parametrized in a proper way. Then we considered 
the equations of motion in the tails of the object, and checked under what conditions the 
linearized equations in the tails can be satisfied. We found, that Klein Gordon PDEs do not 
allow in general for moving breathers (as they do not for stationary ones), only shock waves 
or kinks are allowed. Discretizing these equations, we found that all restrictions are gone, 
and moving breather tails can be generated for any parameters. It is surprising that even 
the nonresonance conditions known to exist for stationary breathers do not appear here. If 
we consider a chain with an acoustic spectrum, then the nonresonance condition reappears 
in some sense, but the forbidden frequency gap shrinks to zero as the velocity is increased 
above the speed of sound. 

Another result is, that we can generate shock wave tails in all cases, with restrictions in 
some cases (acoustic chains) on the velocity. 

IV. NUMERICAL METHODS 
A. The general idea 

Let us consider a chain with the equations of motion given by u n = —dH/du n . Then 
the ansatz (|25|) yields equations of the type (cf. ( p6| ) ) 
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V 2^jM _ 2ikQbV ^M = Q 2 Mz) + f n ({A k ,(z)}, {A k ,{z + n)}, {A k ,{z - n)}) . 

(31) 

The essential feature is that these coupled differential equations contain advanced and re- 
tarded terms. These terms arise due to the interaction on the lattice. Instead of directly 
trying to solve these equations, we consider a lattice governed by the equations 

V 2 A kn (t) - 2ikQ b vA kn (t) = n 2 A kn (t) + £ fn'({A k ' n (t)}, {A k , in+n ,(t)}, {A*, n _ n »(*)}) . 

n' 

(32) 

Here n is again the lattice site label, and with each lattice site n we have an associated 
infinite set of variables {A kn }, k = 0, ±1, ±2, .... Equations (|32|) define a phase space flow in 
the phase space of all variables A kn , A kn . In general trajectories generated by that dynamics 
are not related to solutions of (|3T|). However all fixed points of the map RG t =\ (G t is the 
evolution operator defined by fl3~2|) and R the translation operator that shifts all lattice 
indices by 1) are solutions of (|5T|) . The main reason for that is that all delay and advance 
intervals are integers. 

Once a fixed point (solution) is found, it can be continued using generalized Newton 
methods or steepest descent methods. 



B. Example: DNLS 

Let us investigate numerically moving breathers in the discrete nonlinear Schroedinger 
model (DNLS) with /i = 2 fl2~0|) which is a nonintegrable model. 

The DNLS system has an integrable counterpart which is the Ablowitz-Ladik model 
(ALM) having the following Hamiltonian 

4> n = i(t> n - 1 + t> n+1 )(l + \t> n f) . (33) 

The ALM has moving breathers of the form: 
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_ smh(/i) jkn jn h t ( oa\ 

n ~ cosh(/i(n-W)) ' 1 j 

Here /c and /i are free parameters. 

All other parameters are expressed through k and fi 

= 2 cosh(/i) cos(fc) , (35) 
2 

V = — sinh(/i) sin(fc) . (36) 
A 4 

Here V is the velocity of a breather. 

Now we can test the relations obtained in the previous section. Indeed, using the solution 
([34]) in the tails it can be checked that the relation between the frequency, velocity and 
localization length (p2|),(p3|) is fulfilled. Note that in this case we have only one Fourier 
component which is a specific property of models of the nonlinear Schroedinger type which 
do not generate higher Fourier harmonics with respect to time if a single Fourier component 
is excited. 

Now let us consider the following model, which allows for a continuous tuning between 
ALM and DNLS: 



= i 



^ n _i + # n+1 + |^ n | 2 [(l-a)(^ n _i + ^„ +1 ) + a^g . (37) 



Here < a < 1. 

Let us look for a solution in the form 



= e tUbt g n (t) . (38) 

Then we have 

g n = i g n -! + g n+ i - 2Q b g n + |#„| 2 [(1 - a){g n -i + g n +i) + ocg n ] . (39) 
For a moving breather we have 

g n {t)=g{n-Vt) . (40) 
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The moving breather is a fixed point of the map RGt, T = 1/V. We look for a zero 
minimum of the functional F = \RGtX — X\ where X is a point in the phase space of (p9|). 

We proceed in the following way. First the parameter a is put to zero. The initial point 
in the phase space is chosen to be the ALM moving breather (|34|). Then a is incremented by 
a small value Aa. A minimization of the functional F is performed. Again a is incremented 
by a small value Aa, etc. The algorithm stops when a reaches 1. 

We were able to generate moving DNLS breathers with the value of the minimized 
functional less then 10~ 6 . An example is shown in Fig. 5 (/i = 0.5, k = 1, V = 0.0364). 

Let us mention an important point found during the numerical calculations which is the 
existence of a large (infinite?) number of very (infinitesimally?) close local nonzero minima 
of the functional F near the true fixed point. When a step Aa is made, the algorithm 
minimizes the functional F, and we hope that the value of the functional in the minimum will 
be zero. Surprisingly (and in contrast to the numerical calculations on stationary breathers) 
we found that the value of the functional at the minimumi for a given a is not zero no matter 
how small Aa is. 

On the other hand when Aa is decreased (thereby increasing the computation time) the 
minimum value of the functional at a given a tends to zero, although not being zero for 
any finite Aa. From this we conclude that the true fixed point trajectory is surrounded 
by a dense set of other nonzero minima of the F. The structure of the phase space near 
the moving breather is therefore highly nontrivial in contrast to stationary breathers where 
none of above effects were found. 

These findings are maybe connected to the fact, that the spectrum of Floquet multi- 
pliers of a moving breather has unusual properties as compared to the Floquet spectrum 
of stationary breathers. Namely, the Floquet multipliers of the linearized map around a 
moving breather fill the unit circle densily. Especially there exist Floquet multipliers with 
value +1, which would in general make continuation impossible for stationary breathers. 
The existence of these multipliers can be simply explained. Linearizing the map around a 
moving breather fixed point, we obtain an infinite set of eigenvalues with spatially extended 
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eigenvectors. At large distances from the breather center these eigenvectors will correspond 
to normal phonons. A phonon is given by 

e iQ q t-qn _ ^-q 

It is always possible to cast it into the form 

e m b t e -i(n-vt) ( 42 ) 

with arbitrary numbers f^, V by solving 

Q q = Q b + Vq . (43) 

Indeed we can always find g-values which will do the job. For the case of a stationary 
breather V = 0, and we essentially recover the nonresonance condition, which can be fulfilled 
by choosing Qf, to be outside the phonon band. 

V. CONCLUSION 

In this section we will briefly discuss related work. 

A. Moving kinks 

Moving kinks can be considered to some extend as moving breathers with zero frequency. 
What matters here is that these objects can be represented by one (zero frequency) Fourier 
component in (|3]). Proofs of existence of moving kinks in FPU chains have been obtained 
in |15j by finding them as minimisers of a variational problem and in [plj by analytical 



continuation from the integrable Toda lattice. 



Numerical solutions for moving kinks have been obtained e.g. in [17] and Fourier 
transformations in space are used in the first one, while the second paper treats space- 
periodic solutions, but does not uses Fourier transformations. 
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B. Moving breathers 



There is a large amount of work reporting on moving breathers in FPU chains (e.g. 
|19|). A study of their connection with stationary breathers was started in P0|| . In [^1|] this 
connection was used to numerically obtain moving breathers by exciting pinning modes of 
stationary breathers. 

Finally moving breathers have been obtained numerically for DNLS chains in P^j . 
So far we are not aware of existence proofs for moving breathers. 
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FIGURE CAPTIONS 



Fig.l. Phase diagram for necessary conditions on the existence of moving localized ob- 
jects in Klein-Gordon models in 1 + 1 dimensions (see text). 

Fig. 2 Same as Fig.l, but for the DNLS chain. 

Fig. 3 Same as Fig.l but for a Klein-Gordon chain. The labels M2,M4 indicate the mul- 
tiplicity of solutions in the given part of the parameter space. 

Fig. 4 Same as Fig.l but for FPU chains. 

Fig. 5 A moving breather solution for DNLS (see text). 
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